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Abstract 

The large ratio of color-suppressed tree amplitude to color-allowed one in B ^ tttt 
decays is difficult to understand within the Standard Model, which is known as the 
"B — > TTTT puzzle" . The two tree diagrams contain the up- and charm-quark compo- 
nent of penguin amplitude, Puc, which cannot be separated by measuring B tttt 
decays alone. We show that the measurements of the branching ratio and direct CP 
asymmetry of B~^ — > K^K^ decay enable one to disentangle the Puc with two-fold 
ambiguity. One of the two degenerate solutions of the Puc can solve the B tttt puz- 
zle by giving \C/T\ ~ 0.3 which is consistent with the expectation in the Standard 
Model. We also show that the two solutions can be discriminated by the measurement 
of the indirect CP-asymmetry of B^ K^K^. We point out that if the B vrvr 
puzzle is solved in this way, the corresponding puzzle m B ^ irK decays should have 
a different origin. 
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1 Introduction 



In the Standard Model (SM), rare charmless nonleptonic B ^ nn decays provide valuable 
information on the inner angles of the unitarity triangle of Cabbibo-Kobayashi-Maskawa 
(CKM) matrix [1] . In these decays, dominated by 6 ^ dqq {q — u, d) at the quark level, 
one can measure the angle a through the isospin analysis [2]. In addition, new physics 
(NP) beyond the SM can affect the decay processes [3]. 

In the diagrammatic approach [4], the B ^ ttt: decay amplitudes are dominated by 
color-allowed (suppressed) tree diagram T(C) and QCD-penguin diagram P. These dia- 
grammatic amplitudes are expected to be hierarchical in size, i.e., \C\ ~ \P\ much less 
than \T\. However, the experimental data of the branching ratio (BR) oi B ^ 7r°7r° decay 
requires a large ratio of the color-suppressed tree diagram to the color-allowed one, which 
is in contradiction to the hierarchy. This is known as the "S — > TTTT-puzzle" . 

Recently another S-decay modes, S+ — > K'^K'^ and S° — > K'^K'^, which have identical 
h ^ d transitions at the quark level have been observed by the BaBar [5] and Belle [6] 
Collaborations. They have also measured the time-dependent CP asymmetries, Af and Sf 
defined by 

We study the implications of these newly discovered channels on the extraction of the SM 
amplitudes. Other studies on B ^ KK can be found in [7]. 

We take i) SU (3) flavor symmetry of strong interactions and ii) smallness of annihilation 
and exchange topologies as working assumptions. Then B^ — > K'^K^ and B^ K^K^ 
decays become pure penguin modes. 

The penguin amplitude can be decomposed into two parts depending on the quarks 
running inside the loop, up- and charm-quark penguin P^c and top- and charm-quark pen- 
guin Ptc- In S ^ TTTT decays Puc can always be absorbed into T and C. As a consequence, 

it is impossible to extract pure T or C amplitudes from B — * tttt measurements alone. 

We show that when combining the S — > tttt and B KK decays, we can extract the 
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tree amplitudes and P„c separately. So we can test the ratio \C /T\ without the contami- 
nation of Puc penguin in the b ^ qq transitions. 

We will demonstrate that a solution obtained from the B ^ nn data, B(-B+ — > K'^K^), 
and Acp{B^ K+l{^) has \C/T\ ~ 0.3 which is in accord with the SM hierarchy. First, we 

use the analytic expressions for the BR's and CP asymmetries to show that all the relevant 
amplitudes, T, C, Ptc, Puc, and the angle 7 of the unitarity triangle can be extracted with 
two-fold ambiguity. 

Then we get the diagrammatic amplitudes numerically by performing a x^-fit to the 
current data. Discarding the indirect CP asymmetry of -B ^ K'^K^, Scp{B K^K^), 
which has huge errors, we obtain two degenerate solutions which minimize x^- One solution, 
as mentioned above, is acceptable. However, the other solution gives \C/T\ ~ 1.3 which is 
unsatisfactory. 

Prom the two solutions we can predict Scp{B K^K^). The favored solution with 
\C/T\ ~ 0.3 gives large positive Scp{B K^K^) while the other solution which looks 
unphysical predicts large negative Scp{B K^K^). Therefore if the Scp{B K^K^) is 
measured more accurately in the near future, we can discriminate the two solutions and 
confirm the physical solution. 

It is also well known that \C'/T'\ ~ 1 in 5 — > ttK decays, which is again difficult 
to understand in the SM framework [9]. In B ^ ttK decays, the contamination of P^^ 
does not play an important role because it is suppressed by CKM factors. Therefore the 
"B — > ttK puzzle" , if it remains in the future experimental data, should be solved by other 
mechanisms such as new physics effects. 

This paper is organized as follows. In section 2, using the current experimental 5 — > tttt 
data only, we find that the S — > tttt puzzle is still there. In section 3, we confirm the 
possibility that the B ^ tttt puzzle is simply the effect of Puc and there is actually no 
S — > TTTT puzzle in the "bare" ratio of \C/T\. We conclude in section 4. 
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2 The current status oi B ^ tttt puzzle 

The diagrammatic amplitudes provide a useful parametrization for nonleptonic i?-meson 
decay processes because it is independent of theoretical models for the calculation of 
hadronic matrix elements [4]. The decay amplitudes of S — > tttt's are written as 

V2^(B+ ^ 7r+7r°) = -{T+C + Pew + Pew), 

A(S°^7r+7r-) = -(r + P + lpg^ + E + PA^ , 

(2) 



v^^(S°^7rV) = -(c-P + PEw + lp^w-E-PAy 



Here T, C, P, P^, E, and PA represent the tree, color-suppressed tree, QCD-penguin, 
(color-suppressed) electroweak-penguin, exchange, and penguin-annihilation amplitudes, 
respectively. In general, the diagrammatic amplitudes may have a weak phase as well as a 
strong phase. 

We can further decompose the QCD penguin diagram P which depends on the quarks 
running inside the loop, as 

p = v^dv:,Pu+Vaiv:,Pc+vtdV,iPt 

= Vu,V:,iP^ - + VuV*{Pt - Pe) 

= P„ee'^ + Piee-'^ (3) 

where we use the unitarity property of CKM matrix and explicitly write the weak phase 
dependence for the amplitudes. The relative sizes of the diagrammatic amplitudes are 
expected to have the following hierarchical structure: 
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(4) 



where A is expected to be of order of 0.2 - 0.3. We neglect terms of order A^ or higher and 
thus, consider only T, C, and P. 
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Mode 


B[io-6] 


^CP 


Sep 


B+ 




5.7 ± 0.4 


0.04 ±0.05 




5° 


— > TT'^TT" 


5.16 ±0.22 


0.38 ±0.07 


-0.61 ± 0.08 


5° 




1.31 ±0.21 


0.36 ±0.33 




B+ 
S° 




i.36i°:i 


-0.58lo;gg±0.04 (Belle) 
0.40±0.41±0.06 (BaBar) 


1 Qo+O.SO+O.ll 
-•-■^0-0.73-0.16 



Table 1: The current experimental data for CP averaged branching ratios (BR), direct 
CP asymmetries (^cp) and indirect CP asymmetries {Sep) of 5 — >^ tttt and B —>■ ttK 
decays [5, 6, 10]. 

Taking into account the decay amplitudes containing terms up to 0(A), we find that 
the amplitudes in S — > tttt decays are given by 

-V2yl(B+ ^ 7r+7r°) = {T + C)e'\ 

-x/2^(50 ^ ttV) = C'-' - Ptce-'f - Puce'\ (5) 

where the weak-phase dependence is explicitly written. Prom the above expression we notice 
that the P„c term can always be absorbed into tree diagrams T and C by the redefinition 

f = T±P„e, 

C = C-Puc. (6) 
Then the S — > tttt decay amphtudes are further simplified: 

-- (f±c)e^ 

= fe'^ + Ptce-"", 

-- Ce'^-Ptce-'f. (7) 

As a consequence, the P„c term cannot be detected if we use the P — > tttt data only. 

The current experimental data for the CP-averaged branching ratio (B), direct CP- 
asymmetry (.4cp), and indirect CP-asymmetry {Sep) are shown in Table 1. If we take the 
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-A{B° - 




-V2A{B'' - 
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7 {\T\,Sf) 


i\C\,Ss) 


\Ptc\ 


65.5 ±11.1 (22. 1 ± 0.7. 36.0 ± 22.6) 


(15.8 ± 1.1. -16.5 ± 12.7) 


8.33 ± 1.34 



Table 2: The results for the fit to the B ^ nn data alone. We obtained Xmin/d-O-f = 0.81/1. 
The strong phase of Ptc is set to zero. The magnitudes and angles are in the unit of e^'s 
and degrees, respectively. 

weak phase /3 = 21.2° measured through the b —>■ ccs transition, we have 6 parameters in 
(7). We fit 6 parameters in (7) to the 7 measurements. The results are shown in Table 2. 
The Xmin/d-o.f = 0.81/1 is very good. The angle 7 = (65.5 ± 14.1)° (all the angles in this 
paper are in the unit of degrees) is consistent with the global CKM fit value 7 = 55.6l2!5 
or GS.O^y^Q obtained by CKMfitter group [11]. However, we obtain 

iS = 0.71 ±0.18, (8) 

I I 

which is larger than the expectation (4) by about factor 3. This large ratio is mainly due 
to the unexpectedly large B{B —>■ 7r°7r°). We confirm the B —>■ ttt: puzzle [12] is still there. 



3 Combining B ^ nn with B KK 

Since from the B —>■ ttt: data alone we cannot disentangle the P„c; we include the B —>■ KK 
data in the analysis. There are three decay modes for these b — > d{u)ss transitions: 

A{B+^K+K^) = P + A-^P^^, 
A{B'^K'K^) = P + PA-^P^^, 

A{B^ K+K-) = -E- PA. (9) 

When we neglect small contributions from annihilation and electroweak penguins, we simply 
get 

A{B+ K+K^) = ^(5° ^ K^K^) = P = Ptce"'^(l + re^^-e^(^+^)), (10) 
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where r = \Puc/ Ptc\ and Sue is the strong phase of Puc (We set 5tc = 0.). The A{B^ 
K^K~) is neghgible. Therefore we expect i) B{B^ — > K^K~) is very small, which is 
confirmed by the experimental data [5, 6], ii) B{B+ K+K^) ^ B{B° K^K^), which 
agrees with the data, and in) Acp{B+ K+l{^) Acp{B^ ^ K°'K^). For the last 
relation, the results of BaBar and Belle are contradictory with each other and we need 
to wait for more data to confirm or reject it. As one can see in Table 1, the indirect CP 
asymmetry Scp{B^ K^K^) has still large experimental errors. For these reasons we use 
only B{B+ B{B^ K^K^), and Acp{B+ K+K^) for the numerical analysis 

which will be presented below. 

First, let us demonstrate that the inclusion of just two more data, B{B^ K'^K^) 
and Acp{B^ K~^K^), is enough to extract P^c- Given B{B^ —>■ K~^K^), one can obtain 
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R=— p ' ' ^ = l + r' + 2rcos(/? + 7)cos5„c (11) 



for the ratio of average amplitude squared to |PtcP- The information on \Ptc\ is obtained 
from B TTTT. 

The direct and indirect CP asymmetries for B~^^^^ —>■ K~^^^^K^ are given by 



1 + + 2r cos(/9 + 7) cos 6uc ' 
<c ^rO^^o^n sin 2(/j + 7) + 2r sm{(3 + 7) cos Sue 
Scp{B ^KK) ^ l + r2 + 2rcos(/? + 7)cos<5„e ' ^ ^ 

respectively. When the angles (3 and 7 are given, R and the CP asymmetries are functions of 
r and 5uc only. From (11) and (12), we can determine r and Sue with \Ptc\ and 7 determined 
from B ^ nn decays. 

In Figure 1, we show filled contour plots of ^cp(-B"'"^°^ X+^^^X^) (Figure 1(a)) and 
»Scp(-B° — > K^K'^) (Figure 1(b)). We also draw contours for various values of R with thick 
lines on each plot. For these plots the CP phases are taken to be /3 = 21.2° and 7 = 65.5°. 

In Figure 1(a), the lines for Acp{B+^'^^ K+^^^W) = 0.12 and R = 1.79 which 
correspond to the experimental central values cross at two points denoted by the black and 
white dots. Thereby we get \Puc\ and Sue with two- fold ambiguity once \Ptc\ is obtained 
from the S — > tttt data. 



(a) 



(b) 



Figure 1: Contour plots for constant (a) direct and (b) indirect CP asymmetries as a function 
of r = \Puc/ Ptc\ and the strong phase 5uc (Filled contours). Also contours for constant R 
are shown on each plot. The lines for ^cp(-B^*-°-' K^'^'^^K'^) = 0.12 and R = 1.79 which 
correspond to the experimental central values cross at two points in (a), denoted by black 
and white dots. The points are also shown in (b). 

In Figure 1(b), we can see that the two solutions obtained from Figure 1(a) predict 
very distinct Scp{B^ K^K^). Therefore we can discriminate the two solutions by 
measuring 5cp(-B° K^K'^) in near future. From the figure we note that the prediction 
for 5cp(-B° —>■ K^K'^), especially its sign, is almost insensitive to the change of r unless r 
is too small. 

Now, we proceed to obtain numerical solutions. For this purpose we perform a x^-fit to 
the seven B ^ nn data, i3(5+(°) ^ K+'^^^K^), and AcpiB+ K+W). We do not include 
the CP asymmetries of B^ K^K^ to the fit because they have huge errors and/or BaBar 
and Belle data are not consistent with each other. 

The results for the fit are given in Table 3. Again the quality of fit is very good, 
Xmin/d-O-f = 1.5/2. The weak phase 7 = (65.5 ± 14.1)° is consistent with the global fit 
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1 






65.5 ± 14.1 


(28.7 ±4.8, 25.3 ± 11.6) 


(8.14 ±3.19, -26.12 ±30.6) 


65.5 ± 14.1 


(16.7 ± 1.7, 25.3 ± 11.6) 


(22.3 ±7.1, -9.25 ±9.50) 




( 1 Puc \ 1 ^uc) 


\C/T\ 


(8.33 ±4.33,0) 


(7.89 ±6.82, 174 ± 10) 


0.28 ±0.15 


(8.33 ±4.34,0) 


(6.93 ±3.58, 7.42 ± 11.2) 


1.34 ±0.50 



Table 3: The results for the fit to the B ^ nn and B KK data. We obtained Xmin = 
0.81. The strong phase of Ptc is set to 0. The magnitudes and angles are in the unit of 
eF's and degrees, respectively. 

value by the CKMfitter group [11]. For R, we get R = 1.79 ± 2.07. 

In Table 3, the solutions in the first (second) row corresponds to the black (white) 
solution in Figure 1. For the first solution, the ratio |C/T|, given by 

|C/r| = 0.28 ±0.15, (13) 

is reduced by about factor 3 compared to (8). This agrees with the conventional hierarchy 
in (4). However, as we discussed in section 2, the ratio \C/T\ should be around 0.7. The 
large difference between the ratios \C/T\ and \C/T\ can be explained by the constructive 
interference between C and P„c- As we show in the first solution in Table 3, C and P„c 
are almost same in sizes, while the difference between strong phases Sc and 5uc is around 
TT. Hence the first solution results in the large enhancement in C = C — Puc- The large C 
can accommodate the B{B^ 7r°7r°). In Figure 1, we see that this solution predicts large 
positive «Scp(-B° — > K^K'^). Numerically we predict 

Scp{B^^K^K^) = 0.99 ±0.02. (14) 

Small errors in (14) is because the large errors appearing in r are canceled in the numerator 
and denominator. 

The second solution which gives too large ratio \C/T\ — 1.34 ±0.50, looks unphysical, if 
we take the central value seriously. This can be tested by measuring the <Scp(-B° — > K'^K'^) 
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more precisely. The Figure 1(b) shows that this solution predicts large negative Scp{B^ — >■ 
K^K^). The numerical prediction is 

Scp{B'^ ^ K^K^) = -0.97 ±0.19. (15) 

In [8], the authors could extract the diagrammatic amplitudes from B{B^ K^K^) and 
a theoretical input which is free of the endpoint infrared divergences in QCD factorization. 
The prediction 

S{B^ K^K^) = 0.97 ± 0.02 (16) 
is in full accord with the prediction (14), reinforing the conventional SM hierarchy (4). 



4 Conclusions 

We have considered charmless nonleptonic decay modes, — > tttt and B KK. We find 
that the large ratio of color-suppressed tree to the color- allowed tree which is known as 
the — > TTTT puzzle" is mainly due to the up- and charm-quark penguin P„c contributions 
which are always included in the tree diagrams in the form of T = T -|- P^c and C = 
C — Puc- The Puc can be determined with two-fold ambiguity by using the measured 
B{B+ K+K^) and ^cp(^^ ^ K+K^). The two-fold ambiguity can be hfted when 
<Scp(-B° — > K^K'^) is measured with more accuracy. Consequently, the bare ratio \C/T\ 
without the contamination of Puc can be determined. 

Conversely, we can predict Scp{B^ K^K^) by imposing the conventional hierarchy 
to the solutions. For the theoretically preferred value \C /T\ ~ 0.3, we get 

Scp{B^ K^K^) = 0.99 ± 0.02. (17) 

If the B ^ 7171 puzzle disappears in this way, the B 7iK puzzle which is represented 
by \C'/T'\ = 1.6 ± 0.3 [13] will become more prominent. Its solution will require other 
mechanisms in the SM or new physics beyond the SM. 
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